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^ ■ Abstract 

Qh| We calculate radiative corrections to the Casimir effect for the massive complex scalar 

field with the A(/> 4 self- interaction in d + 1 dimensions. We consider the field submitted 

to four types of boundary conditions on two parallel planes, namely: (i) quasi-periodic 

boundary conditions, which interpolates continuously periodic and anti-periodic ones, (ii) 

Dirichlet conditions on both planes, (iii) Neumann conditions on both planes and (iv) 

mixed conditions, that is, Dirichlet on one plane and Neumann on the other one. 
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^§ 1 Introduction 

> : 

In 1948, H. B. G. Casimir P] computed for the first time the energy shift of the vacuum state 

of QED caused by the presence of two parallel perfectly conducting plates close to each other. 

43 ■ 

For simplicity, and as a first approximation, Casimir considered non-interacting fields 1 and that 

the conducting plates influence only the (quantized) electromagnetic field. For this case, the 

vacuum energy shift can be obtained by the shift in the zero-point energy of the electromagnetic 

field. 

Zero-point energy is not a peculiarity of the electromagnetic field, so that any relativistic 

quantum field shall exhibit a Casimir effect when submitted to boundary conditions (BC). A 
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1 The expression "non-interacting field" along this paper shall always mean that the field does not interact 
with itself or with any other field. However, it is sensitive to the presence of material plates, but the interaction 
between the field and the plates shall be simulated by idealized boundary conditions (BC). 



detailed discussion about many aspects of the Casimir effect can be found in ref's El HI El IEj 
and references therein. 

However, any meaningful quantum field theory describing nature is constructed by inter- 
acting fields. Therefore, in order to study the Casimir effect in its totality we should take into 
account all orders in perturbation theory. In terms of Feynmann diagrams, we should compute 
the influence of the BC's not only in the one-loop vacuum bubbles, but also in the two-loop 
ones, and so on. 

Although there is a vast literature concerning one-loop contributions to the Casimir effect 2 , 
the same is not true, comparatively speaking, for two-loop contributions. In the context of 
QED, two-loop calculations can be found [TJ] and for scalar fields in [To ^ ITd^ IT71 IT%1 IT31I2U] . 

Though the first radiative corrections (two-loop contribution) to the QED Casimir effect 
does not have an experimental motivation, since it is a\ c /a smaller than the one-loop effect, 
where a = 1/137, A c is the Compton wavelength of the electron and a is a typical distance in 
Casimir experiments, it may be relevant in the QCD bag model |13j . as pointed out in jE]. 

However, there is indeed an important theoretical motivation for studying the Casimir 
effect at two-loop level, namely. To test renormalizability of QFT under the influence of BC 
at higher orders in perturbative theory. Idealized BC at higher orders may lead to unsolvable 
inconsistences. 

In what concerns interacting fields, one of the most popular models of Quantum Field Theory 
in the literature is the so called A0 4 self-interaction model. This occurs due to its simplicity 
and its applicability to describe many physical phenomena, as spontaneous symmetry breaking 
and Bose-Einstein condensates. 

Although the Casimir effect for the scalar field with A0 4 self-interaction has been studied 
recently, there are still a few questions about this subject we want to consider here. The first 
one is the fact that, for some situations involving Neumann BC's, the first radiative correction 
to the Casimir effect diverges for d = 2 spatial dimensions in the massive case [T%1 IT§1 l2"U] . The 
second question, related to the first one, consists in the necessity of introducing renormalization 
surface counter-terms for the scalar field in some situations to get finite results, as pointed out 
in the work of reference ^Hj • 

This paper is devoted to the study of the first radiative corrections to the Casimir energy 



2 The most popular method for computing the one-loop Casimir effect was that introduced by Casimir in 
1948 I], based on the zero point energy of the field. However, many other equivalent techniques are available for 
this calculation, as for example Green's function method :.6;, generalized zeta- function method [7j, Schwinger's 
method 001113111] and other methods (EJ. 



for the massive complex scalar field 0, with the well known self-interaction A0 4 . Specifically, 
using euclidean coordinates, (xq, Xi, ..., xa —'■ z), in D — d+ 1 dimensions, and a system of units 
where % = c = 1, we can write the Lagrangian density of the field in the form 

£E = |<V>| 2 + m 2 |0| 2 + A|0| 4 + £ c t (1) 

where C ct is the lagrangian counter-terms. We consider the field always submitted to boundary 
conditions on two parallel planes, located, in our coordinates system, at Xd = z = and 
Xd = z = a. 

In fact, the A0 4 theory is renormalizable only for d = 2, 3 spatial dimensions, even when 
the field is not submitted to any boundary condition. Therefore, our calculations will be valid 
only in these cases (d = 2, 3). 

This paper has the following structure: in section 2 we develop, up to order A, the general 
quantum theory for the complex massive scalar field described by the Lagrangian density (|T|) 
and submitted to BC's at two parallel planes (or lines in d = 2 spatial dimensions) located 
at Xd = z = and Xd = z = a. In section 3 we study the case of quasi-periodic BC, which 
interpolates continuously the periodic and anti-periodic conditions on one spatial direction. In 
section 4 we study the case where the field obeys Dirichlet BC at both planes. In section 5 we 
consider the case where the field obeys Neumann BC at both planes, and finally in the section 
6 we study the case where the field satisfies Dirichlet BC at Xd = z = and Neumann BC at 
Xd = z = a. The last section is devoted to some conclusions and final remarks. 

2 General theory 

Despite the fact that we are considering a field submitted to BC's, its Quantum Theory follows 
in strict analogy to the case where there is not any BC's. The energy of the vacuum state of 
the field can be written as [T§1 |2*T] 

rT/2 

-T/2 

where x = (x±, ..., Xd — z), xq is the temporal coordinate and Ce is the Lagrangian density (JTJ). 
It is worth emphasizing that we are using Euclidian coordinates. 

For all boundary conditions we considered in this paper, the renormalization of the quantum 
theory described by the Lagrangian (PQ) is assured, in order A, by the Lagrangian counter-terms 

C ct = 5m 2 \4>\ 2 - ( Cl 5(z) + c 2 6(z - a)) |0| 2 + 6A , (3) 
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E = lim In 



! V(f)V(j)* exp ( - / dx f d d x£ E (x , x) 



(2) 



where <5A is a variation in the cosmological constant due to the A0 4 interaction (i.e., the change 
in the vacuum energy which is due solely to the interaction, and not to the confinement), 5m 2 
is the mass renormalization counter-term, and c\ and C2 are surface renormalization counter- 
terms. These last terms can be interpreted as extra renormalizations for the field mass, which 
must be taken only at the planes Xd = z = and Xd = z = a. 

By means of standard perturbative methods of quantum field theory, equation (J2J) gives the 
first radiative correction for the Casimir energy per unit of area, 



e 1 , r a 

= g 1 = dz 

A Jo 



2\Gl(x, x) + 5m 2 G c (x, x) — [ci5(z) + C2<5(z — a)]G c (x, x) + 5A 



(4) 



with E 1 representing the first radiative correction for the Casimir energy, A = j dx±... J dxd-i 
is the planes area and G c (x,y) is the Green's function of the field submitted to the boundary 
conditions, but in the absence of the A0 4 self-interaction term. In writing expression (0J) we 
discarded divergent contributions that come from the integration in the exterior regions of the 
planes. These contributions do not depend on the distance a, and then do not contribute to 
the Casimir force. 

In order to calculate the integrals presented in (@J), we have first to write the Green's function 
G c (x, y) in its spectral representation, which can be done by considering the field modes between 
the planes. Since all boundary conditions considered here are imposed on the Xd = z coordinate, 
these field modes can be written as 

(f)(x) = $(x±)(p n (z) , (5) 

where we defined x± = (xq, x%, ..., Xd-i), &(x±) has the usual expansion for free fields and <p n (z) 
satisfies the BC's. In addition, the functions <f n {z) form an orthonormal set according to 

( Pn(z)(p* nl (z) = 5 n , n , , (6) 



with 5 ntn i designating the Kronecker delta. 

Now, the spectral representation for G c (x,y) reads 

cm =iS?^ ^ ■ ^ - *>] ? $$$& p> 

where k± = (k , ..., kd-i) and q n is the <i-th Fourier factor, quantized by the boundary conditions. 
By analytic continuation, equation (J7|) gives 122] 

G c (x, x) = j^T (l - fj £ dt\ n (z)&(z) , (8) 



with 

^n = \/ql + m 2 . (9) 

Although the summation written in (jHJ) is divergent in the general case, the expression (JSJ) 
shall allow us to calculate the integrals presented in (@J) by analytic continuation in all situations 
considered in this paper. 

The renormalization counter-terms are determinated from the Green's function of the field 
(with boundary conditions) which, in first order in A, reads 

G(x, y) = G c (x, y)-J d D £G c (x, £)£ C (£)G C (£, y) , (10) 

where 

S c (0 = 4AG c (e, + Sm 2 - [ Cl 8fa) + c 2 8(£ d - a)] (11) 

is the self-energy of the field in the presence of the BC's. 

Although SA does not depend on the boundary condition imposed on the field, it is conve- 
nient to calculate this counter-term for each case separately. This is done by imposing 5A to 
be a-independent and setting to zero the contribution for the first radiative correction of the 
Casimir energy which is proportional to the spatial volume between the planes, that is, the 
contribution linear in the distance a and in the area A. 

The mass counter-term 5m 2 is calculated by imposing that the quantity E c (£) satisfies the 
following conditions: 

^c(0 == ^ finite between the planes 
a lim£ c (0 = 0, (12) 

except, possibly, on the planes. In addition, we consider 5m 2 to be a-independent. 
From the image method, we can write G c (x, x') in the form [T9~|l27]] 

G c (£,O=G (£,O+Gc(£,O (13) 

where Go(£, £') is the Green's function for the non-interacting field without BC's, and (j c (£, £') 
is a correction induced by the BC's which is finite in the region between the planes, but can, 
possibly, diverge on the planes. Besides, G c (£, £') has the property vanishing when the distance 
between the planes goes to infinity. 

Substituting (p"3|) in fTTf we can see that we are taken to a divergent term proportional to 
Gq(£jO- I n or der to have accordance with the first condition (|12|) . this divergence must be 



canceled. It is done by taking 

4AG (£,0 + W = //, (14) 

where /z 2 is a real finite constant. Substituting (JHJ) into (fTTj) we have 

E c (0 = 4AG C (£, + fi 2 - [aSfa) + c 2 5(U - a)] . (15) 

Taking the limit a — > oo of the above expression and considering the second condition (fT2|) we 
have that /i 2 = (where we used that lim a ^ 00 G c (x, y) = 0). Therefore, equation (JT3J) yields. 

where we used the Furrier representation of Go(£, £) in euclidean coordinates and calculated 
the above integral by analytic continuation [22] • 

The mass counter-term expressed in (llfjjl is valid for all boundary conditions and dimensions 
considered in this paper. 

The surface counter-terms c\ and c 2 are determinated by imposing that the first order 

correction to the Green's function G(x, y) written in (jlOJI is finite between the planes. This is 

equivalent of imposing 

f d£ d G c (x, 0S c (0G c (e, V) =► finite . (17) 

Jo 

We can not analyze the above expression in a general case. It must be considered for each 
specific situation separately. 

Now, let us calculate the first radiative correction to the Casimir energy for each BC's 
mentioned previously. 

3 Quasi-periodic boundary condition 

In this section we consider the case where the field satisfies the so called quasi-periodic boundary 
condition, namely, 

0(x_i_, Xd = z = a) = exp(i9)(j)(x±, Xd = z = 0) (18) 

where 9 is a parameter that can vary continuously from to n. Note that for 6 = and 9 = it 
we have, respectively, periodic and anti-periodic conditions. 

With conditions (|T5|) the functions (f n (z) written in (jSJ) take the form 

1 ( 2n7r + ° \ 

(Pn{z) = —=exp[i z\ , n = 0, ±1,±2, ... , (19) 

'a \ a I 



with the factors q n appearing in (J2J) given by 



Qn 



2nir + 6 



a 



(20) 



For the BC's (|18|l there is no need to consider the surface counter-terms c\ and c 2 in the 
Lagrangian density 0, because they are not helpful to renormalize any divergence at all 3 . 
Therefore, equation (JIJ) becomes 



>(i) 



(.(,./, , i 



2AGg(x, x) + 5m 2 G e (x, x) + <5A 



(21) 



where Gq(x, x') designates the Green's function for the non- interacting <fr field submitted to the 
condition (fTHj) . 

Using equations (jHJ), (jUJ), ([HI) and (|2T?j) we can write 

r(i-d/2) 



Gf)(x,x) 



(47r) d / 2 a n 



E ^ (d < 1), 



with 



uv, 



f2nn + 9 s 



\{- 



+ m 2 . 



(22) 



(23) 



Inserting expression ()22j) into equation (|2*T|) and arranging terms, we obtain 



S, 



(i) 



2A 



T(l-d/2) ~ d „ 2 a^m 2 
~~ n 4A 



+ a 



5A 



(5m 



2N2 



JA 



(24) 



( 47r )d/2 

In order to compute the summation that appears in equation (|24jl it is convenient to reex- 
press it as 



" d-2__ f^L) d ' 2 v (tlA ™ l_ 

a 



E - 

n=— oo 



2 ' 2vr ' 2tt ' 



(25) 



where we used equation (|23|) and defined the function X> by the expression 

/ n \ 2" 



p hs R£ 



V+ [ n+ 2i. 



Re(s) > 1/2. 



(26) 



This function has the analytic continuation to the whole complex s-plane given by 



2 T' l/ '2* 



\ \/7ri/ 



l-2s 



r( s ) 



I ! -s - - I + 4 ^ cos {n6) -f . 

£=i (n7rz/) 1 /2 . 



(27) 



3 Recall that imposing BC means, in some since, to compactify one spatial dimension. Consequently, we 
are left to a compact manifold without boundary and hence it is meaningless to talk about surface terms. 



which exhibits simple poles at s = 1/2, —1/2, —3/2, . . . Inserting equations (J25J) and (|2*7|) into 
equation (|21jl we have 

2A f am d - 1 



£ 



(i) 



a \(4vr)( rf + 1 V 2 
(5m 2 ) 2 



r [ !___- ] + 4 V cosM) ^-dM™™) 
^ 2 J + ^ l j (nma/2)(^)/ 2 



+ 



a £m 2 1 

1a - | 



+a 



5A 



2 3 A 



(28) 



In order to fix 5 A we require that linear term in a in the above expression of £} must vanish, 
as explained before. So we have 



5A 



{5m 



2\2 



2 3 A 



(29) 



Inserting (|T6|) and (|2"9"|) into equation (J28|) we finally arrive at the first radiative correction 
to the Casimir energy per unity area for quasi-periodic boundary conditions: 



sPia) 



Xm 



d-l 



2d-2 7T d+l a d-2 



oo 

y^ cos {n6 

n=l 



K(d-i)/2(man) 

' n (d-l)/2 



(30) 



As mentioned before, expression (|3()|) is meaningful only for d = 3, 2 spatial dimensions; 



-(I) 



£?\a) 



Xm 2 



™ K^nma) 
2 , cos(rw) 

n=l 



n 



siV 



(31) 
(32) 



5=3 27T 4 a 

— — *-- [lnfl + e" 2ma - 2e" ma cos(0))l 2 
^=2 a (4?r) 2 L V v >)\ 

where we used expression ()113|) . 

For the special cases of periodic {9 = 0) and anti-periodic {9 = n) BC's, equations (|3Uj) 
reduce to 



•(i)/ 



-(I) 



£^{a) = £^ (a) 



Xm 



d-l 



•(i)/ 



-(I) 



2d-2 7T d+l a d-2 

Am d_1 



EVlfl) = £&(*) = p,-,^,,-. 



E 

ra=l 

oo 

E 



K {d _ 1)/2 {man) 

n (d-l)/2 



(33) 



% 1)/2 (2man) 
~^ (man)'''" 1 '/ 2 



i+l)/2 y Vl)/2M 



n=l 



(man)( d 1 )/ 2 



Similarly, equations (J3TJ1 and (|S2*|) take the form 



£?{a) 


£^{a) 


£ ( p\a) 


4%) 



d=3 



d=3 



d=2 



d=2 



^■0=0 \ a , 



(1) 

9 = 7 

(1) 
9=C 

£il(«: 



*firto 



£&(a 



Xm 2 

d=3 ~ 2vr 4 a 

Xm 2 



=3 2vr 4 a 
A 



rf=2 2vr 2 a 
A 



^=2 2vr 2 a 

8 



E 

.71=1 

OO 

E 

_ra=l 

[ln(l- 

ln(l + 



K\ {man) 



n 



Ki{2man) — Ki{man) 



n 



(34) 

(35) 

(36) 

(37) 
(38) 



If we further take the limit m — > in equations (JHDJ), (jS)) and (|32|) we have 

2 



& ( V 



sPia) 



sPia) 



A 



d-1 



m=0 



d=3 

m=0 



d=2 
m=0 



27r d+l a 2d-3 



A 



2vr 4 a 3 

A 

2^ 



E 

n=l 



cos nt 



??, 



E 

vn=l 

_A_ 

2^ 



cos(n#) 



/■?' 



d-l 



Ms. 



In 2 sin (0/2 



(39) 
(40) 
(41) 



where we used equations (jll6|) and (J117J) . and B 2 (x) = x 2 — x + 1/6 is the Bernoulli polynomial 
of second degree [21j . 

For the specific cases of periodic [6 = 0) and anti-periodic [6 = n) conditions, equation (J3~9"|) 
leads to. 



S?(a) 



£%\a) 



m=0 



m=0 



4%(a) 



Sipja) 



A 



'd-l 



m=0 



m=0 



27T d+l a 2d-3 

A 

27r d+l a 2d-3 [ ~^pi 



1 



The above results for d = 3 become 



S?(a) 



£%\a) 



d=3 
m=0 



d=3 
m=0 



«,(<*) 



'(9= 

£Er(a 



\t\d-l), 



i 



d=3 
m=0 



d=3 
m=0 



2 3 3 2 a 3 ' 
1 



which could be obtained from expression (|4T)|) . Similarly, for d = 2, we get 



4 1J (a) 



si'V 



d=2 
m=0 



d=2 
m=0 



*£>(«) 



dlfo^ 



d=2 
m=0 



d=2 
m=0 



OG 



A 'H2f 



2n 2 a 



(42) 
(43) 

(44) 
(45) 

(46) 
(47) 



which could also be calculated using expression (|4~T|). 

The above results are in perfect agreement with those presented in the literature J7]. It is 
worthwhile to emphasize that for m = and d = 2 the first radiative correction to the Casimir 
energy per unity area diverges in the periodic case (8 = 0), as expressed in equation (J4*6|) . A 
similar situation occurs for the field at finite temperature: there, for d = 2 spatial dimensions, 
we have a divergence when m = 0. 

From result (J3*0"|) we can see that the first radiative correction to the Casimir energy per 
unity area for quasi-periodic boundary conditions is always positive, whatever parameter 6 we 
use. 
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4 Dirichlet boundary conditions - DD 

In this section we consider the field submitted to Dirichlet boundary conditions at both planes, 
it is 

4>(x±, Xd = z = 0) = (f>(x±, Xd = z = a) = , (48) 

a situation which, from now on, we shall refer to as DD conditions. In this case, the functions 
ip n appearing in (jHJ) are given by 



and the q n factors become 

As a consequence, equation (J£]J) gives 



\2 fnnz 

- sin 

a V a 



n=l,2, 



mr 

q n = — 

a 



U n 



n7T Y , 2 

a ) 



(49) 



(50) 



(51) 



Once the field satisfies DD conditions, the corresponding Green's function Gdd shall, like- 
wise, satisfy (J48J1 . So, their contributions to the integral in (jlj) vanish, whatever values of Ci and 
c 2 we use. Therefore the counter-terms c\ and c 2 can be set to zero and equation (JJJ) becomes 



4d= / dZ d \2\G 2 DD {£,Z) + 8m 2 G DD {ti,£) + 8h 

Jo L 



(52) 



where G\dd (£,£') is the Green's function for the non- interacting field submitted to the DD 
conditions. Besides, taking into account the Lagrangian counter-terms (jHJ), we can see that the 
contributions of the counter-terms c\ and c 2 for the action (integral of ((IJ) is always zero, due 
to the condition (|4*K|h 

Using equations (|HJ) and (J49)L we can write the Casimir energy ()52j) in the form 



f 1 

°DD 



2XT 2 (l-d/2) 
a (47r) d 

+ 5m 



5> 

L \n=l 



rf-2 



1 oo 

+^E^ 4 

Z n=l 



2 r(i- W 2) )v - y . 



(to)* 2 



oo 

E-: 

n=l 



(53) 



Completing the square, we have 



f 1 

°DD 



2A 
a 



T(l - d/2) ,52, „_„ <5mV 2 r2 



( 47r )rf/2 



E -n" 2 + 



71=1 



4A 



r 2 (i - d/2) » 2d -A 

2(4vr)^ £< » 
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2\2* 



+ a \8A - 



(5m 
2 3 A 



(54) 



By the same arguments presented for the quasi-periodic BC's, the counter term 5A for 
DD conditions is also given by (f2U|) . The summations appearing in (|H4"|) can be calculated by 
analytic continuation (see the Appendix). Therefore, using equations (fTHjl . (J2EI), ((HI} , ((EH) and 
(jllOj) . we have the first radiative correction for the Casimir energy per unity area with the field 
submitted to DD BC's in d spatial dimensions: 



S 1 



DD 



2A 
a 



Aam° 



lL V K (d-V/2( 2man ) 

( 47r )(d+l)/2 A, ( man )(^l)/2 



+ 



+ 



4am d " 1 r(l-d/2) 

2 (4 7 r)( rf + 1 )/ 2 2(4vr) d / 2 r " 

2am 2d-3 p2( 1 _ d / 2 ) 



,'-2 y> K (d ^ 1)/2 (2man) 



E 



man)'^ 1 )/ 2 
-^(2d-3)/2(2man) 



( 47r )d+i/2 r(2-d) ^ (man)( M ~ 3 )/ 2 
where we discarded an a-independent term which does not contribute to the Casimir force 
For d = 3, 2 spatial dimensions, we have, respectively 

Am 2 

- | J - 



(55) 






DD 



d=3 



°DD 



32% 2 a 

A 7 2 



2 ~ K x (2man) 

n=l 



-1 



/? 



(56) 



( 1 + - ln(l - exp(-2ma)) J 



d=2 32vr 2 a 

where we used the fact that T(— 1/2) = —l^pix and T(— 1) 
result (JH7|) we used equation (|114|l and the limit 

T(l-d/2) 



1 



(57) 



ex) for the result (JEEJ)- For the 



lim 

d^2 



r(l-d/2 



r(2-d) 

with 7 being the Euler's constant. 

For a vanishing mass, the limit m — > in equation (JH5J) yields 



27 



(5f 






DD 



,3-2d 



m=0 



2A 



( 4 ^)d+l/2 



d- 1 



+ 



Wl^V 



(4vr) 1 /2 ^ 2 
^2d-3\C(2t/-3) 



C 2 (d-1) + 



(59) 



, 2 y v 2 / r(2-d) 

where we considered that lim m ^om d ~ 2 = 0, used equation (jllfjjl and £ designates the Riemann 

zeta function. Taking the limits d — > 3, 2 of equation (|59|) . we have, respectively 

A 

.;. -:-,,„ ...() 2 7 3 2 a 3 ' 

a r 12 

"'--*■ "'-" 2 5 3a vr 



C D_D 



6 DD 



"(271 + 7 2 ) 



(60) 
(61) 
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with 



2~£hf 



7l=n Q „ 2 



8=1 



(62) 



[(l-s)C(s)\ 

Results (jnOj) and (JHT|) agree with those presented in the literature for m = (THl IT7j . 
Expression (jfjQjl could be obtained from (JHfij) in the limit m — > and with the aid of (|118j) . The 
limit m — > of equation (JH7|) is divergent. For this case (m = and d — 2) we have to take 
first m — > considering lim m ^ md ~ 2 = 0, and after that we take d — > 2, as we have made for 
computing (jfiTj) . 



5 Neumann boundary conditions - NN 



For this case, the field satisfy the following conditions 



d(j)(x) 



dx 



d(j)(x) 



x d =z=0 



dx 



. 



(63) 



Xd=z=a 



With the restriction (|53|l the functions ip n appearing in (jSJ) read 

<fn(z) -- 



'2-Snjo (nnz 
cos - 



n = 0, 1,2,3,..., 



V a \ a I 

with the parameter q n and frequencies o; n , defined respectively in (J2J) and (0), given by 



(64) 



n7r 
a 



a 



(65) 



For NN conditions the first radiative correction for the Casimir energy per unity area (|2} 
becomes 

4iv = f <*& [2AG^(e, + <WGW(£, - c[<*(&) + <*(& " «)]GW(£, + <5Al , (66) 

JO l J 

where we used the fact that c\ = c 2 due to the symmetry of the system, and defined c := c\ — c 2 . 
Using equations (jSJ) and (|51jl. expression (jo^jl takes the for 






2AT 2 (l-(d/2)) 



NN 



(4tt) 



1 CO 

E^ 2 ) + ^ 2 

\n=0 / Z n=l 



+ 



+ 5m 



a r(i-(d/2)) 



E ^ + 



cT(l-(d/2)) 



a (47r) d / 2 ;f 



£(2-5 n , )^n" 2 + ^A 



(67) 
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Completing square and recalling that uq = m, what can be seen from (|fj5j) . we write equation 
()67j) in the form 



f 1 

°NN 



2A 
a 



T(l-(d/2)) 



( 4vr )rf/2 



E ^ + 



5m 2 a 



n=0 



4A 



+ 



+ 



T 2 (l - d/2) 



2(4tt) 



5> 

n=l 



2d-4 



+ T 



c T (1 - d/2 ) / 
A 2(47r) d / 2 



m d " 2 -2E^- 2 

\ n=0 / 



+ 



+ aUA- 



(5m 



2\2' 



2 3 A 



(68) 



In order to calculate the surface counter-term c, we consider condition (fTTj) and write the 
self-energy £jvat(£) expressed in ffTTjl for the NN conditions as, 



Z NN (0 = 4AGW(£, + <5™ 2 - c[S{U) + <K& ~ a)} . 



(69) 



In this case, integral (fT7j) diverges due to a divergent contribution which comes from Eatat(£) 
and is given by the Green's function evaluated in the coincident point limit and also lying at 
one of the planes, that is, Gnn(£,, with £ d = z — or £ d = o. This divergence of the integral 
(|I7jl coming from £aw(0> expressed in (jHEl), can not be renormalized by the mass counter-term 
8m 2 . 

Therefore, the divergence of the integral (|17)l comes from the integration over regions near 
the planes, and condition (|T7j) is equivalent of imposing 

d£dG NN (x,Z,)Y lNN (Z,)G NN (Z,,y)+ d£, d G NN (x,£)E NN (£)G NN (£,,y) -> finite, (70) 

for an arbitrary e, since £ < a/2, and for i,?//^. 
Using the fact that 

gg^fcg = o , (71) 

we see that Gnn(x,£,) does not depend on the coordinate ^ near the planes. Therefore, for a 
small e, equation (|70|l can be written as 



G NN (x, l €d= o [ / d£ d E NN (£, f )] GW(£, 1/) | €d= 



+ GW0r,Ole d = a [/ dfdEjvjvfoOjGWfo^l^a -> /mite . 
Once Gnn(x, £)L =0a is finite (for ^ / i 11 = z), we have that 

/ d£dE NN {£,£)+ d£ d E NN (€,€)^> finite. 

JO Ja-e 



(72) 



(73) 
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Adding the finite contribution J e a e T,nn(£,£) to both sides of (J73|) . we conclude that condition 
(fl7J) . for the case at hand, can be substituted by he following one: 



5W&0- finite. (74) 

Substituting (J69|) into (f?4"J) and using expressions (jHJ), (jHljl and (J65|) we obtain 



r(l-(d/2) °° . 
4A V, .,.„ ! V w„ + Wa - c -> /mzie , (75) 

(4tt)<V 2 ^ 

where we considered f£ 5(z d ) = J a 5(z d — a) = 1/2. 

Taking into account the definition of u n , expressed in (|53jl . the analytic extension (J110|) and 
the result (|16|). condition (|75|) can be written in the form 



rfl-(d/2)) , 



+ 16A am d ~ l T K V-W&™™) _» finife (76) 



Due to the presence of the T function, the first term on the right hand side of (J75J) produces 
a divergent a-independent contribution for d = 2. This problem can be circumvented by setting 

c = 2A F \~ n y{ 2 ' m d - 2 - Xm d - 2 a NN (d) , (77) 

where (iNN^d) a dimensionless function of the dimension d and the product m d ~ 2 <TNN(d) is 
finite for any d. 

In contrast, for d = 3, expression (|7fi)) does not contain any divergence and the use of the 
surface counter-term c is not necessary. Therefore, awN^d) introduced in equation (|77jl is chosen 
such that c = for d = 3: 



ri-(3/2) i 

However, this is not sufficient to fix <JNN{d) completely. Hence, for d = 2, the value a(d = 2) 
can not be determined. 

The last counter-term appearing in (jHHj) is the cosmological constant counter-term SA. It 
is given by (J2HJ) and determined by the same arguments used for the boundary conditions 
considered previously. 

Substituting (fTHjl . (J2HJ) and (J77j) in (jHHjl . using the definitions (jHHjl and the analytic extensions 
(JllOj) and (jllOj) . discarding an a-independent term which does not contribute to the Casimir 
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force and collecting terms, we have the first radiative correction to the Casimir energy per unity 
area for NN conditions, namely, 



£ l 



NN 



2A 
a 



^am d - 1 ~ K {d „ 1)/2 (2man) 
( 47r )(d+i)/2 Z, ( mcm )(d-i)/2 



+ 



4am d - 1 /T(l-d/2) <r(d) N 



+ 



( 47r )(d+i)/2 ^ 2(47r) d / 2 2 

2am 2d-3 r 2^ _ d / 2 ) ~ K {2d ^ )/2 (2man) 



d-2 y K {d - 1)/2 (2man) | 

n=l 



fman)^- 1 )/ 2 



(4vr) rf + 1 /2 r(2-d) 



E 

ra=l 



(man) ( M ~ 3 )/ 2 



(79) 



Taking d = 2 and d = 3, which are the relevant cases, we obtain 



°NN 



f 1 



NN 



d=3 



d=2 



Am 2 
32n 2 a 

w 

32n 2 a 



' 2™ K x {2man) Y ' 



7T 



n=\ 



n 



JO) 



4ttct(2) 

7 



j +-ln(l-exp(-2ma)) ~( 1 + 



4ttct(2) 

7 



oo and used equations (|75j). (J55|) and 



where we considered that T(— 1/2) = —2^/lx, V{— 1) 

dm. 

The zero mass limits are obtained by taking m — > in equation (|79|1. The calculations are 
performed analogously to those made for the DD case (jo^j) and the result is given by 

„3-2d 



f 1 

°NN 



m=0 



f 1 
°DD 



m=0 



2A- 



l7r )d+l/2 



d 



(47T)V2 



|C 2 (d-i)+ 



+ r 2 



'2-d\ (2d-3\((2d-3) 



52) 



For d = 2, 3 we have 



f 1 

°NN 



A 



f 1 



AW 



d=3,m=0 



d=2,m=0 



2 7 3 2 a 3 ' 

A 
_ 253^ 



12/ 

1-^(271+7' 



(83) 
(84) 

and (l84l) are 



where 71 is defined in (162)1 and 7, as usual, the Euler's constant. The results 
in perfec agreement with those found in the literature [T6|ll7j. 

The result (183)1 could be obtained by taking the limit m — > in (|80j) . However, if we take 
this same limit in ()81)) we get a divergent result, in contrast with the finite result given by ( |84jl . 
The reasons for such a disagreement are the same as those presented for the DD case, where 
an identical situation occurs. 
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6 Mixed boundary conditions - DN 



In this section we consider the field satisfying Dirichlet condition at z = and Neumann 
condition at z = a: 



(fi(x ± ,x d = z = 0) 



dcj)(x) 



dx 







x c i=z=a 

which restrict the tp n functions written in (J5J) to be given by 



¥n{z) 



sm 



n + 



1\ nz 



n = 0,l,2,3,... . 



Also, q n and u; n , defined respectively in (J2J) and (jHJ, take the form 



9n 



1\ 7T 

n+- - 
2/ a 



U-'„ 



n+ 2)a. 



m z 



55) 



(86) 



57) 



By the same reasons as those used for the DD conditions, the counter-term c\ appearing in 
expression (J1J must be zero. As a consequence, the first correction to the Casimir energy (J1J 
is written in the form 



Or 



DN 



di d 2\G 2 DN (£, + 6m 2 G DN (£, - a$(& - a)G DN (£, f) + 6A 



where, for simplicity, we have done c = c^- 
Using equations (jSJ) and (|8*5jl. the integral 

2Ar 2 (l-(d/2)) 



yields 



°DN 



(47r) c 



1 



Vn=0 / Z n=0 



2d-4 



+ 



+ 5m 



a r (i-(d/2)) 

(47r) ^ ;i u 



E^" 2 + 



c r(i-(d/2)) 



E ^ + «^ 



9) 



a (4tt)<*/ 2 /( „ 

The surface counter-term c is calculated considering the self-energy £ expressed in (jlljl . as 
well as the condition (|17|). In this case, these expressions are given, respectively, by 

Sdtv(0 = 4AGW(f , + Sm 2 - c5(£ d - a) , (90) 

/ di d G DN {x,Cj^DN{0 G DN{y,Cj ^ f inite ■ ( 91 ) 

Jo 

As in the case of NN conditions, the integral (|91|) is divergent due to the integration over 
the regions near the planes, £ = 0, o, so condition (|9*T|l can be substituted by the equivalent one 

/ d£ d G DN (x,£)T lDN (£,£)G DN (y,£)+ d£ d G DN (x,£)T, DN (£,£)G DN (y,£) -^ finite , (92) 

Ja-e JO 
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where e is an infinitesimal positive arbitrary parameter. 

The second term on the left hand side of (JQ2} is finite. This can be shown by considering 
that it is a good approximation to take ^ ~ in this term, once the integration runs over 
< t;d < e and e is infinitesimal. But, when £ d ~ 0, it can be shown that Gdn(x,0 ~ (,d and 
also T*dn{£.) ~ (£d) 1-d - Therefore the second term on the left hand side of (f9*2"|) reads 



f dz d G DN (x, ox DN (OG DN ( y , o ~ f di d u^y-% ~ f d& (&) 

JO jo jo 



~< / ^£d — finite , (93) 

Jo 

where, in the second line, we used that d = 2, 4. Equation (|93j) states that the second term on 
the left hand side of (}9~2*j) is finite, and therefore, condition (J9~2"|) reads 

f d&GWfo 0S DW (e, Z)G DN (y, - /imte. (94) 

Ja—e 

Using the fact that 

. (95) 



dG DN {x,£) 



Zd=a 



we see that Gdn(x,£) does not depend on £<$ near the plane ^ = a. Therefore, for a small e, 
condition (jHljl reads 

G D n{x,£)\u=o. [J_ d^DN^,0) G D N(y,0kd=a -> finite . (96) 

Taking into account the fact that (jdjv(£, Okd=a is finite, we have 

f d£ d E DN (Z,Z)^ finite. (97) 

Ja—e 

Adding the finite contribution f^7 2 e d^^DN '(£,£) and introducing a convenient well behaved 
function (the cosine), the condition (|9Tj) reads 

/■a I 7TZ \ 

/ dz d T, DN (z, z) cos — > finite . (98) 

Ja/2 \ a J 

The introduction of the cosine in (J98)) is justified by the fact that, along the integration 
domain, it is a finite function with constant (negative) sign, and has the property of not 
depending on ^ for £ d ~ a. Hence, it does not disturb the divergences which appear in (|9*7j) . 

Substituting (|9T)|) in (|9*H|) and using equations (jHJ), (|HH|) and (jHZJ) we have 

r(i-(d/2) 



4A- 



(4^)^/2 



1 1 °° 

! -1-2 , - 1 V^ , ,d-2 



r£< 



4^0 +-2^ w - 

^ " n=0 
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5m 2 — > finite , (99) 

7T 2 



where we used that /°/ 2 d£d$(£d~ o) — ~ 1/2- Taking into account the definition of u n expressed 
in (J%7|) . the analytic extension (jlllj) and the result (fTHjl. condition (JOH)) reads 



-A 



r(i-(d/2) 



•16A am 



d-l 



(4vr)(d+i)/2 



2E 



( 4vr )rf/2 

K(d-i)/2(4man 



2a 



(d-2)/2 



C 

+ 2 



+ 



E 



•K"(d_i)/2(2man) 



/irate . (100) 



^i (2man)( d - 1 )/ 2 ^ (mem)^- 1 )/ 2 

Analogously to the NN case, the T function in the first term on the left hand side of (|1()()|) 
produces a divergent a-independent contribution for d = 2, which is canceled out by taking 

r(l - (d/2)] 



2A- 



m d - 2 - Xm d - 2 a DN (d) 



101 



(47r) d / 2 

where a^N(d) is an arbitrary dimensionless function of d, the product m d ~ 2 ar,N(d) is finite and 
we took into account that c can not depend on the distance a. 

For d = 3 the T function in (|100|) is finite and the introduction of the surface counter-term 
c is not necessary. Then we must take 

r(i - (3/2)) = _ j_ 

2vr 



cr(d = 3) 



(102) 



(47r) 3 / 2 

what makes c equal to zero in (jlOlj) for d = 3. 

Similarly to the NN case, for d = 2 the function <JDN{d) can not be determined. 

The SA counter-term present in (fKHj) is given by (f2HJ) and it is determined by the same 
arguments used in the previous sections. 

From equations (fTHjl. (f2TJ|) . (J%7|) . (J%HJ) and (jlOlj) . and also, using the analytic extension (jlllj) . 
we obtain 



°DN 



2A 
a 



4am 



d-1 



2E 



AT(d-i)/2(4man) 



E 



iC( d _i)/ 2 (2man) 



+ 



( 47r )(d+i)/2 ^ ^ (2 ma „)(rf-i)/2 ^ (man)^" 1 )/ 2 

ra-(d/2)) ,_,/ra 



+ 



2A 

a 



4(4vr) d / 2 

T(l- 



(d/2)) md _ 2 



(d/2)) md _ 2 



(d/2)) \ 

r(2 - d) ) 

T(l - (d/2)) 



771 



2 cr DN (d) 



+ 



771 



2 c^Z)Af(d) , 



(103) 



4(47r) d / 2 "" V r(2-d) ~; I 

which is the first radiative correction to the Casimir energy in d dimensions for DN conditions. 
For d = 2, 3 dimensions we have, respectively, 

Am 2 



f 1 

°DN 



d=3 



32vr 2 a 



2 ~ Kidman) - K 1 {2man) 



7T 



n=l 



n 



(104) 
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f 1 



DN 



d=1 



A 7 2 

327r 2 a 



2 (l — exp(— 2ma)\ 

— In ; 

7 Vl-exp( 



Zj//fclXJ I 

{—2ma) J 



1 + 



27RX gJV (2) \ 



+ 



/ | 27rgggr(2) y 



(105) 



where we used that T(— 1/2) = —l^fii and T(— 1) ^ooas well as equation (J102)) for the result 
IjITHj) . For the result (fTUHjl we used equations (JSHJ and (fTTlj) . 
The zero mass limitof equation (J103J) is given by 

'd- 



£ 



DV 



m=0 



2Aa 



3-2d 



(^ 



2-d 



_(4tt)C«H-i) 

where we used equation (|llfi|) and the fact that lim r , 
For d = 2,3, equation ()l()fi|) yields 

£dJV 

&DN 



l)T z 



,m 



\ 2 J 

d-2 = 0> 



C 2 (^-i) 



(106) 



A 



d=3,m=0 



d=3,m=0 



2 9 3 2 a 3 
A 

2 3 7r 2 a 



In 2 (2) 



(107) 
(108) 



where for the result ()107j) we used that T(l) = 1 and £(2) = vr 2 /6, and for the result (|108|) we 



used that 



lim 

d^2 



(2 2 -° 



1 r 5 



d 



C 2 (d-1) 



47T 



(109) 



(47r)(« f + 1 ) v ^ " V 2 

Results (|l()7jl and (J108J) agree with those found in the literature [TH] . 

Similarly to what happened to the previous cases, result (J107J) could be obtained by taking 
the limit m — > in (|l()4j) with the aid of ()118jl . The same procedure does not work for d = 2 
once this limit when applied to equation ()105|) diverges. For d = 2 we must take first m — > 
and then d = 2 as we did to obtain (J108J) . 

7 Conclusions and Final Remarks 

In this paper, using dimensional regularization, we calculated the first radiative correction to 
the Casimir energy of a complex massive scalar field with a self-interaction term A0 4 in 2 + 1 
and 3 + 1 dimensions. The restriction of these specific dimensions is justifyed by the fact that 
for higher spatial dimensions the A0 4 theory is not renormalizable, even without boundary 
conditions. 

We considered the field satisfying boundary conditions at two parallel planes. Specificaly, we 
studied four tipes of boundary conditions for the field, namely: quasi-periodic conditions (which 
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interpolate continuously the periodic and anti-periodic cases), Dirichlet conditions, Neumann 
conditions and mixed conditions. In the zero-mass limit we obtained the results previously 
published in the literature |IBJ HZl CHJ- It is interesting to note that for a massless field the 
results for DD and NN conditions are equal for both cases of 3 + 1 and 2 + 1 dimensions. 

This was not expected since any two-loop contribution to the Casimir energy depends on 
the form of the field modes and not only on the corresponding eigenfrequencies as it occurs with 
the one- loop contribution to the Casimir effect, which can be obtained through the zero-point 
energy of the field. 

However, for the case of a massive field, the two-loop contribution to the Casimir effect are 
different for DD and NN conditions, as it can be seen from equations ([56)1 . (J57J1 . (jHUj) and (|8ip. 

For configurations involving Neumann conditions in 2 + 1 dimensions, one needs surface 
counter-terms, and the first radiative correction to the Casimir energy can not be completely 
determined. There remains an undetermined arbitrary factor introduced by the surface counter- 
terms. Maybe this fact is an indication of the limitation of idealized boundary conditions (in 
this case, Neumann condition) commonly used in the literature, or even a limitation of the 
meaning of 2 + 1 models. 

The appearance of the undetermined function of the dimension <jnn(cI) alerts us for other 
kind of problems that may appear when computing two-loop diagrams under Neumann bound- 
ary conditions. Maybe the regularization/renormalization program can not be achieved success- 
fully with such an idealized boundary condition. This problem deserves further investigation 
and we think it is not a peculiarity of self-interacting scalar fields. It seams that first radia- 
tive corrections to the QED Casimir effect with infinitely permeable plates instead of perfectly 
conducting ones has unremovable divergences 4 . 
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Appendix 

In this appendix we calculate summations needed to obtain some results presented in the paper. 

4 This problem is under investigation at the moment, but preliminary results suggest that the two- loop 
Casimir energy of QED with permeable plates is ill defined 27 
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It is well established in the literature that the analytic continuation of the so called Epstein 
function is given by [2EJ 



E 

ra=0 



m 2 4 I — I 



a J 



-s/2 



-m 
2 

1 -s 

-m + 



tE 



n=— oo 

am l ~ 



nil 



nr + | — J 
's- 1 



2-i -s/2 



+ 4£ 






2" ' 20rT(s/ 2 

which is valid for any s complex, except s = 1, —1, —3, —5, ..., where it has simple poles 
From the previous equation, we also have the following analytical extension 



(110) 
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ra=0 



-. \ 2 / \ 2-i -s/2 
1 \ / 7T X 



m +n+ 2 U, 



E 



n=0 



m 2 + — 
\2o. 



-s/2 



E 

n=0 



nn 



m 



-s/2 



am 



rfcl + 



+ 



2y/jfT{s/2) V 2 
2am 1 ~ s 



v^r( s /2) 



~ Ay s)/2 (4man) 
£i (2man)( 1 - s )/ 2 



-E 



if(i- s )/ 2 (2raan) 



~^ (man)' 1 " 8 '/ 2 
Another important summation used in the paper is calculated using the fact that 



:nr 



2tt 
K 1/2 (2man) = l-^^ e ~ 2man 
l V Iman 



(112) 



which allows us to write 

Ki/2(2man) 



DO 

y^ cos(n6 
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y 

(man) 1 / 2 Ama ~^ 



p —n(2ma—i0) „—n(2ma+i0) 



+ 



n 



n 



TV 



E 



da(2ma — i9)e 



-cm(2ma—i9) 



+ / da(2ma + iO)e 



-an(2ma+i9) 



Ama^ 

/yj- /-OO °° /"OO °° 

" (2ma-i6) daY e- an{2ma ~ ie) + {2ma + iO) da V 

„— a(2ma- «#) 



Ama 



Ama 



{2ma — iff) I da 



, 1 p— a(2ma— i9) 



f°° e 

(2ma + 16) / da 

Jx 1 — 



-an(2ma+id) 



a(2rna+i9) 



p—a(2ma-\-i9) 



Ama 



— lnf 1 + e~ 4ma - 2e" 2ma cos(fl) 

7/7. V 



In the special case where 9 = we have 

~ K 1/2 (2mcm) 
~^ (man) 1 / 2 



2ma 
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7/7. V 



-2ma 



;ii3) 



(114) 



Along the text, in the zero mass limits, we used the result J2E] 



hmi^)-»— ^ , (i/>0) (115) 



in order to write 



lim 



2v 

m 



K u (2man)~\ x '"'--' 



oo 



^2T(u)vi 



n=l 



2v f x 

{2manY (man) 1 ' 



2 v /,v 7 a 

where £ is the Riemann zeta-function. Analogously, we have 

fC;, (man 



ir(z/)c(2i/)i , (lie) 



lim 



oo 

m^ J^ cos(n 



71, 

71=1 " 



^rM^E^. di7) 



a "l w 



In the zero mass limit we have 



~ KA2man) 1*1 1 , . ?r 2 , . 

^i man ~~ 2(ma) 2 ^ n 2 " 2(ma) 2 ^ 1 ; 12(ma) 2 ' l ; 

where we used that Ki(x) = (l/x), as x — > 0. 
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